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THE  INITIAL  VALUE  PROBLEM  FOR  Till  EQUATIONS  OF  MOTION 


OF  COMPRESSIBLE  VISCOUS  AND  HEAT-CONDUCTIVE  FLUIDS 
Akitaka  Matsumura*  and  Takaaki  Nishida** 


Si.  Introduction  and  Theorem 

The  motion  of  the  genera]  isotropic  Newtonian  fluids  is  described  by  the  five 
conservation  laws: 


1,  2,  3 


p.  +  (pu  )  «  0 

t  x . 

3 

(1.1)  <  U^  +  U1Ui  +  —  P  =  ~  {  (p  (U*  +  u3  ) )  +  (p  ‘  u3  )  )  , 

it  X.  p  X.  p  X.  X.  X  X.  X, 

1  3  i  3  31 

0P 

u-’e  +  — u3  =  -L  { (<e  >  +  y)  , 

v  c  X.  PC  X.  PC  X  X. 

3  3  3  3 

12  3 

where  p:  density,  u  =  (u  ,  u  ,  u  )  :  velocity,  0:  absolute  temperature, 

P  =  p(p,  0):  pressure,  p  =  p  (p  ,  0):  viscosity  coefficient  u'  =  U '  (p  .  0):  second 
viscosity  coefficient,  <  =  k  (p  ,  0)-.  coefficient  of  heat  conduction,  c  ■  c  Ip  ,  0): 

p ' <u'  )  :  dissipation 


heat  capacity  at  constant  volume  and  'C  =  (u3  +  u^  )  “  +  p1  (u^  )': 

2  x,  x.  x. 

*  3  3 

function.  We  consider  the  initial  value  problem  for  (1.1)  with  the  initial  data 


u.2)  (p ,  u,  0)  (0,  x)  =  (p 0 .  uQ,  0Q)  (x)  ,  x  e  F3  . 

We  seek  the  solutions  in  a  neighbourhood  of  any  constant  state  (p ,  u,  0)  =  (p ,  0,  0), 
where  p,  0  are  any  positive  constants.  Thus  we  assume  the  following  natural 
conditions  on  the  system  (1.1)  of  hyperbolic-parabolic  type  throughout  this  paper: 

(i)  p,  c,  p,  p*  and  k  are  smooth  functions  in  0  =  {  (p ,  u,  0):  ) p  -  P  |  ,  |u|, 

1 0  -  e |  <  t). 

2 

(ii)  3p/3p ,  3p/30  >  0,  c,  p,  k  >  0  and  p1  t  -  P  >  0  in  0, 
where  e  <  min  {p ,  0). 
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i'ii!>t  n>wi  ite  the  system  (1.1)  by  the  change  of  (lie  unknown  and  known  variabl* 

as  follow*  i  p  •»  p  t  p,  u  "  U|  8  **  0  ♦  8,  p(p  ♦  p,  8  +  0)  *♦  pip,  8),  U  (p  *  p«  u,  5  <  0) 

►*  g  (P ,  u ,  8 )  and  so  on . 

s  0  -  1  i 

1.  (p  ,  u)  p  ♦  (p  ♦  p)u  ♦  u  p  “0 

xj  x3 


n.i > 


whore 


(l.t) 


L*  (u)  u*  -  ii  u'  —  (v>  +  M  * )  u  ' 

t  x.X  *.*. 

j  3  i  ) 


L  (0)  8  -  K  8 

t  XjXj 


l  -  r 

G  -p  p  -  P.  0  +  g 

p  x .  8  x . 

i  i 


i  -  1,  2.  3 


-P 


u’  + 


3  x 


i 

l 


g 1  -u'  u*  +  In  lu'  +  u '  )  +  g  ‘  (u'  ) )/ (p  +  p ) 

*J  *}  Xj  Xi  Xi  *1 

g"*  -u'  0  +  <x  8  +f)/(0+  p)c 

x,  x.  x. 

}  1  1 


More  wo  also  use  the  abbreviations 


g  -  g/(p  +  o)  ,  ii’  »  u’/tP  +  P)  ,  pp  ■  p  /(P  1  8 > 


p  -  pft/(P  *  ’  Pj  -  (8  +  0)  Ve/(P  *  (>)c  and 

x  *  x/(p  ♦  P)c  . 

bet  H  (l  »  1,  2 ,  •  •  •  ,  5)  be  the  Sobolev  space  with  the  norm  II  II ;  of  the 

L.,- functions  having  all  tho  t-th  derivatives  of  functions.  Define  i  by  the 

Sobolev's  lemma  so  that  for  II  fll  ,  <  t  wo  have  max  |f|  <_  C  II  fW  ,  <  i  .  Denote 

.  o.  a ,  n, 

D  f  {8  f/3  x,1  3  x,  3  x,  for  all  a,  a  +  a  ♦  a  »  t) ,  !  »  1,  •  •  •  ,  -1. 

12  3  I  «  j 

Tho  initial  data  for  (1.3)  are  given  by 


(t.’j) 


(p,  u,  0)(0)  -  lp(),  u(1,  8^)  <  H  n  Lj  for  <  -  3  or  4 


t 


;  E)  for  some  E  <  r  , 


£ 

Tho  solution  is  sought  in  the  space  of  functions  X  (0,  *• 
v  *•  3  or  4,  where  for  O'-  t  ^  <  t  ^  <  «*> 

x'  (t  t,  lE)  .((p,  u,  0)  (t)  :  0  (t,  x)  t  ^(tj,  t  ,  I  H*)  n  01  (t  ^ ,  t  t  !  H* “ 1 } 
ul(t,  x)  ,  0  It,  x)  t  C°<V  )  HC)  n  Cl(t1,  t.,  i  Ht_2)  n 

(1.6) 

£  +1 

n  W  *2  1  "  ’  •  i  -  1.  2,  3  .  ami 

*3 

sup  it  (p.  u,  o)(t.)ii*  ♦  /  lip  (8)11 ;  +  ii  (u.  exsHi2  ,  ds  -  e‘J  . 

t  <t<t„  t  1 

1 - 2  1 

Theorem 

Consider  the  initial  value  problem  (1.3)  (1.5)  and  let  the  initial  data  have  the 
noi  m  for  t  ■ 


<*•  7>  E  -  II  (p ,  u,  0)  (0)||  +  II  (p ,  u,  O)  (0)11  <  »  . 

£  L1 

Then  there  exist  positive  constants  and  CQ  «  •  (C0  4Q  t)  such  that  if 

Kf  <  50'  then  tho  problem  (1.3)  (1.5)  has  t  he  unique  solution  (p ,  u,  0)  (t)  in 
large  time  such  that 

(p,  u,  0)  (t)  r  X*  (0,  •  i  C0  E{) 


and  it  li.i-.  I  lie  decay  rate 

(1.8)  ||  (p. 


U,  0)  (till  ,  <  C  E./(l  ♦  t)  3  4  . 

a  "*  0  V 


In  particular  if  u ,  p  •  and  K  do  not  depend  on  p .  then  tho  above  assertion  holds 
for  $  «  3  also. 


In  II)  we  obtain  the  same  type  of  result  in  the  more  restricted  case  of  a 
poly tropic  gas.  We  also  refer  the  reader  to  (1)  for  the  bibliography  of  other  known 
results. 


-3- 


5 2 .  Proof  of  Theorem 


estimates  for  the  solution  in  X 
Theorem  2.1  (local  existence) 


Theorem  is  proved  by  a  combination  of  a  local  existence  theorem  and  a  priori 

t 


Consider  the  initial  value  problem  (1.3).  t  >_  t  ^ .  Let  the  initial  data 


(p ,  u,  8)  (tj)  t  H  for  £  “  3  or  4 


Then  there  exist  three  constants  5^  >  0,  <  “  (C^  6  ^  <  e )  and  t  >  0  which  are 


independent  of  t  such  that  if  II  (p ,  u,  0 )  ( t  111  <  6  ,  then  the  problem  (1.3) 

1  1  £  1 


(1.8)  has  the  unique  solution 


(p  ,  u,  8)  (t)  t  X  (tj,  tj  +  T  ;  0^  II  (p  ,  u,  8 )  ( t j ) II  ) 


The  proof  for  £  -  4  is  the  same  as  that  for  poly tropic  gas  in  11).  We  need  an 


approximation  of  the  initial  data  in  H  and  the  L,  energy  estimate  for  the  case 


1-3. 

Theorem  2 ■ 2  (a  priori  estimates) 

Suppose  that  for  the  initial  data  having  the  norm  E.  <  •»  for  £  *  4,  there  is 


a  solution 


(p ,  u,  8)  (t)  £  X  (0,  T  ;  E) 

for  some  T  >  0  and  some  E  <  e.  Then  there  exist  positive  constants  e  (<  c), 


6 ^  and  (C  y  5.,  <  e)  such  that  if  E  <  c ,,  and  E^  <  then  the  solution  has  the 


a  priori  estimates 


(P,  U,  6l(tl  f  X  (0.  T  (  C,  Ej) 


Proof  of  Theorem 


Take  SQ  -  min  U^,  «2,  tj/Cj,  Sj/Cj.  Cj/Cg  C2>  and  CQ  -  C,.  We  use  the 
standard  continuation  argument  of  local  solution  on  10,  n  t ) ,  n  ■  1,  2,  •  •  •  to  get 
the  global  solution.  In  fact  by  the  local  existence  theorem,  the  definition  of 
*0  JnJ  t>ie  dssumFtio"  <Et  *  $Q)  we  have  a  positive  constant  x  and  a  local 
solut ion 

(P,  u,  6)  (t)  c  X1  (0,  !  |  Cj  Ejl  . 

BV  Et  (  Cj  <  r2  and  Eg  <  4Q  <_  4^,  a  priori  estimates  gives 

(P,  u,  6)  (t)  t  Xv  (0,  T  i  C„  Ee )  . 

But  by  Eg  <  4^  4^  and  the  local  existence  theorem,  we  have  again 

(p,  u,  0)  (t)  c  X*  (T,  2t  ;  Cj  C,  Eg)  . 

Now  by  C2  Eg  <  Cj  C2  4Q  c2  and  Eg  <  8Q  4n,  a  priori  estimate  shows 

(p,  u,  6)  (t)  e  X*<0,  2t  ,•  C,  E  )  . 

2  { 

Thus  we  can  continue  the  same  arguments  on  (n  T,  (n  +  1)t)  and  10,  (n  +  1)t) 
successively  n  »  2,  3,  •  •  •  . 


-5- 


<  l.  A  Priori  Estimates 


Lcnnu  i.  1 


(i)  X^  depends  on  i  |c|  only  and  X^  =  0  iff  |t |  “  0,  j  «  1,  •  •  •  ,  4. 

(ii)  X^  ^  Xu ,  j  f  k,  for  all  |d  except  at  most  four  points  of  | £ |  >  0. 

J  K  '  '  - - -  ■  ~  - - 

liil)  There  exist  positive  constants  r^  <  r 2  such  that  X  has  a  Taylor  (Laurent) 
senes  expansion  for  |t|  <  rx  (|c|  >  r  ,  respectively).  Specifically,  the 
Taylor  series  has  the  form 


>/a2  +  b2  i  If!  +  - 


2  .2 


2  - 


a  ■»  b  )  (2p  ■»  p ' )  ■*•  b  k  , . 


2  2 
2  (a  +  b  ) 


(i  Id)  + 


X2  =  X*  (complex  conjugate) 


(3.4)\ 


X,  =  *2.*.  (i  |c|)2  +  a-2  fa2  b2M2w  +  u-)  -  a2  <)  (.  |5|}4  + 

(a2  *  fa2,' 


3  2^2 
a  +  b 


=  li  <i  Id 


(iv)  rank  (X4  -  A)  =  3  for  all  |d  >  0. 


(v)  The  matrix  exponential  has  the  spectral  resolution 

(3.5) 


4  tX.(U 

^  =  l  e  2  P.  (C) 
j=l  3 


(3.6) 


(3.7) 


for  all  |c|  except  at  most  four  points  of  |?|  >  0. 

IIP.  (C)ll  <  C  for  k|  <  r,  . 

]  -  -  1 

It  has  the  estimate  by  the  modification  of  the  right  hand  side  of  (3.5)  near 
the  points  of  a  multiple  eigenvalue 

lle“(5)ll  <  C(1  +  t)3  e_8t 


for  |c|  >  and  a  positive  constant  8. 


-7- 


Ia  mma  ) .  2 


There  is  a  constant  C  •  C(t)  such  tliat 


(3.8) 


ii  f  <  u  >  M  .  iiF(u)ii  <  c  it  mi; 

h  “  2 


ID  F (U) H  <  C  BUS  IHW  for  k  »  1,  2  . 

—  2  k*2 


In  particular  in  the  case  of  (1.9) 


(3.9) 

P^>£0!n^on3j3 


ii d"  f (0)11  <  c  nmi  mull,  ♦  ii u.  en.) 

2  a  9 


There  exrst  6^,  and  such  that  if  E^  <  6^  and  E  <  e^,  then  U(t) 

satisfying  (3.2)  has  the  estimates 


(3.10) 


Him)!  i  c3  Et  (1  +  t) 


-3/4 


\ 

l 


/  II  U (t)ll  ^  ds  <CjE(  , 


where  t  =  4  in  general  and  i  ■  3  for  the  case  (1.9). 

The  Proposition  is  a  consequence  of  Lemmas  3.1  and  3.2.  In  fact  we  have 


II  U  ( t )  II  <_  CQ  E£  (1  +  t) 


-3/4 


C  /  (1  +  t  -  s) 


-3/4 


ii  U  (s) II  2  ds 


V. 


II  D  U  ( t ) II  <  C„  E.  (1  +  t) 


-5/4 


-  0  t 


C  /  (1  +  t  -  s) 


-5/4 


II  U  ( s )  II  . 


•  (IIU(s)ll,  +  II U (s) II  )  ds  ,  k  =  1,  2  . 

2  4 

3/4  7 

Therefore  for  M(t)  •  sup  (1  +  s)  IIU(s)ll  ,  we  have  M(t)  C  E,  +  C  M(t)  + 

0<s<t  1 

C  E  M(t),  where  E  is  the  norm  (1.6)  assumed  on  the  solution.  Thus  we  get  the 

conclusion  of  Proposition  for  f  =  4. 

Next  we  have  to  obtain  the  estimates  for  the  higher  derivatives,  which  is  given 
by 


-8- 


i  ami  C  such  that  if  E  <  6  ,  E  <  t 


L(t)  ;dtut  u?i>  the  estimates  1 1 . 1)  ,  then  On  energy  estimates  hold 


Since  vie  have  already  obtained  the  necessary  estimates  (3.10)  for  the  lower  order 


derivatives  of  the  solution,  the  proof  of  Proposition  3.4  is  given  successively  with 


respect  to  k  and  m  in  the  same  way  as  that  for  polytropic  gases  in  (11.  In  fact  let 


nonlinear  terms  g  in  the  right  hand  side  of  (1.3) 


Lemma  3.5 


We  have  the  estimates  for  k  =  0,  1 


The  estimate  (3.11)  for  k  =  2  is  given  by  the  integration  on  x  t  K 


the  equality 


Integrate  by  parts,  use  (3.10)  and  Lemma  3.5.  The  estimate  (3.12)  for  m  =  3  is 


obtained  by  the  integration  on  xe]R,0f_t^_T  of  the  equality 


Integrate  by  parts,  use  the  equation  (1.3)  and  (3.10),  (3.11)  for  k  =  2  and  Lemma  3.5 


We  can  proceed  to  get  (3.11)  for  k  =  3  by  (3.14)  and  so  on.  The  detailed  arguments 


using  the  Friedrichs  mollifier  and  the  estimates  for  composite  functions  are  the  same 


as  that  in  11 J,  and  so  we  omit  them  here 
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